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ABSTRACT 

The EOLE Experiment is revisited to study turbulent processes in the lower stratosphere circulation 
from a Lagrangian viewpoint and resolve a discrepancy on the slope of the atmospheric energy 
spectrum between the work of Morel and Larcheveque (1974) and recent studies using aircraft 
data. Relative dispersion of balloon pairs is studied by calculating the Finite Scale Lyapunov 
Exponent, an exit time-based technique which is particularly efficient in cases where processes 
with different spatial scales are interfering. Our main result is to reconciliate the EOLE dataset 
with recent studies supporting a k^^^^ energy spectrum in the range 100-1000 km. Our results also 
show exponential separation at smaller scale, with characteristic time of order 1 day, and agree 
with the standard diffusion of about 10^ m^s^^at large scales. A still open question is the origin 
of a k~^^^ spectrum in the mesoscale range, between 100 and 1000 km. 
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I. INTRODUCTION 



The EOLE project consisted in the release of 483 constant- volume pressurized balloons, in 
the Southern Hemisphere mid-latitudes throughout the period September 1971-March 1972, 
at approximately 200 hPa. The life-time of these balloons was from a few days to about one 
year, with an average value of about 120 days. Their motion was basically isopycnal except 
for small diurnal volume variations of the envelop of less than 1% and inertial oscillations of 
a few meters in the vertical, excited by wind shear and small-scale turbulence. The position 
of the balloons and meteorological parameters were periodically transmitted to satellite by 
ARGOS system. The trajectories of the EOLE experiment still provide nowadays the most 
extensive data set of experimental quasi-Lagrangian tracers in the atmosphere for observing 
the properties of the medium-to-large-scale motion at the tropopause level. 

Both Eulerian and Lagrangian analyses have been performed by several authors. Morel 
and Desbois (1974) deduced the mean circulation around 200 hPa from the balloon flights, as 
formed by a mid- latitude zonal stream with a meridional profile characterized by a typical ve- 
locity ~ 30 ms~^ inside the jet, overimposed to meridional velocity field disturbances of much 
smaller intensity, ~ 1 ms""*^, and to residual standing waves acting as spatial perturbations 
of the zonal velocity pattern, producing the typical shape of a meandering jet. These results 
have been largely confirmed by operational analysis since then. Morel and Larcheveque 
(1974), hereafter ML, investigated the synoptic-scale turbulent properties. They measured 
the mean square relative velocity and the relative diffusivity of balloon pairs, and found 
essentially two major regimes for Lagrangian dispersion: exponential separation for time 
delays less than 6 days and standard diffusion for larger times. These authors also observed 
that the scaling of the relative diffusivity with the separation length between two balloons 
agreed with a direct 2D turbulent cascade, with energy spectrum E{k) ~ k~^, or steeper, 
in the range 100-1000 km. Further Eulerian analyses of large-scale velocity spectra by Des- 
bois (1975) were compatible with the scenario proposed by Morel and Larcheveque (1974) 
about isotropic and homogeneous 2D turbulence with a k"^ energy distribution up to scales 
~ 1000 km. 

Later, other authors reached for different conclusions after observing energy spectra in 
the low stratosphere, measured from experimental data recorded from commercial aircraft 
flights. Gage (1979), Lilly (1983), Nastrom and Gage (1985). Their picture suggested a 2D 
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turbulent inverse cascade, characterized by the E(k) ~ k spectrum, inside the interval 
of scales 100-1000 km. 

More recently, Lindborg and Cho (2000, 2001a and 2001b) computed velocity spectra 
using data recorded during the MOZAIC program and also found a k~^^'^ spectrum. They 
suggested a dynamical mechanism different from 2D inverse cascade where energy is injected 
in the large scales by breaking Rossby-gravity waves and generates a chain process down to 
smaller scales. Their hypothesis is supported by the observation of a downscale energy flux, 
whereas 2D inverse cascade should exhibit upscale energy flux (Lindborg and Cho, 2000). 

We wanted to reconsider this issue by performing a new analysis of the relative dispersion 
properties of the EOLE balloons within the framework of dynamical system theory. Relative 
dispersion properties are analyzed through the computation of the Finite-Scale Lyapunov 
Exponent, or shortly FSLE (Aurell et al. 1997, Artale et al. 1997, Boffetta et al. 2000a). 
The FSLE is based on the growth rate statistics of the distance between trajectories at fixed 
scale, and is a better tool at analyzing scale dependent properties than plain dispersion, 
as explained below. This new method has been already exploited for studies of relative 
dispersion in atmospheric and oceanic systems (Lacorata et al. 2001, Joseph and Legras 
2002, Boffetta et al. 2001, LaCasce and Ohlmann 2003, Gioia et al. 2003) and also in 
laboratory convection experiments (Boffetta et al. 2000b). 

This paper is organized as follows: in section El we describe the FSLE methodology; 
section IIIII contains the results obtained from our analysis of the EOLE experimental data; 
in section HVl we discuss the physical information that can be extracted from this paper and 
possible perspectives. 

II. FINITE-SCALE RELATIVE DISPERSION 

Generally speaking, most flows exhibit a range of scales over which fluid motion is ex- 
pected to display different characteristics: a small-scale range where the velocity variations 
can be considered as a smooth function of space variables; a range of intermediate lengths 
corresponding to the coherent structures (and/or spatial fluctuations) present in the veloc- 
ity field over which velocity variations are rough but highly correlated; a large-scale range 
over which spatial correlations have definitely decayed. In each of these ranges, relative 
dispersion between trajectories is governed by a different physical mechanism (chaos, turbu- 
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lence, diffusion) which can be, in principle, identified from the observations. In fully devel- 
oped three-dimensional turbulence, motion is only smooth under the Kolmogorov dissipative 
scale. In the free stratified atmosphere (above the planetary boundary layer), turbulence is 
a relatively rare event: motion is most often smooth but for some localized convective or tur- 
bulent events, associated with mesoscale systems, that mix momentum and tracers. Hence 
one expects to find a smooth (chaotic) dispersion range ended at a scale characteristic of 
the spacing of mixing events, followed by a range covering the mesoscale to synoptic range, 
and finally standard diffusion at planetary scale. This view is supported by the ubiquitous 
observation of long-lived laminated structures in the free troposphere (Newell et ai, 1999). 

In order to fix some terminology, we will use both symbols R and 6 for indicating the 
distance between balloons: the former will be considered as a quantity function of time, the 
latter as an argument for scale-dependent functions. 

A. Diffusion and Chaos 

Diffusion is characterized in terms of diffusion coefficients, related to the elements of a 
diffusion tensor defined as 



where Xi{t) and Xj(t) are the i-th and j-th coordinates at time t, with i,j = 1,2,3. The 
average operation () is meant to be performed over a large number of particles. The diagonal 
elements are just the diffusion coefficients. When the D^s are finite, then the diffusion is 
standard. This means that at long times, after the Lagrangian velocity correlations have 
decayed (Taylor, 1921), the variance of the particle displacement follows the law: 



In presence of a velocity field characterized by coherent structures, it is more useful to 
observe the relative dispersion between the trajectories, rather than the absolute dispersion 
from the initial positions, given by ()II.2|) . which is unavoidably dominated by the mean 
advection. 

In the case of the EOLE experiment, where observing the expansion of balloon clusters 
with more than two elements is a rare event (ML) a measure of relative dispersion is given 



A, = lim -((x,(t) - (x,(0)))(x,(t) - (x,(0)))) 

t — >oo Zt 



(ii.i) 




(11.2) 
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by the mean square inter-particle distance: 

(i?2(t)) = (||x(-)(t)-xW(t)|p) (II.3) 

averaged over all the pairs (x'^™-', x*^")), where m and n label all the available N trajectories. 
Notice that the norm in pi.2|) and pi.3|) must be defined accordingly to the geometry of the 
fiuid domain, i.e. in the atmosphere we use the arc distance on the great circle of the Earth 
connecting the two points. The quantity in pi.3jl can be measured for both initially close 
pairs, balloons released from the same place at short time delay, and so-called chance pairs, 
balloons initially distant which come close to each other at a certain time and, later, spread 
away again (ML). Consistency of the average in eq. pi.3|) requires all the trajectory pairs 
to have nearly the same initial distance, a condition which strongly limits the statistics. At 
long times, {R(tY) defined in eq. (III.3|1 is expected to approach the function ADt, where the 
4 factor accounts for relative diffusion. When it happens that (i?(t)^) ~ t'^'^ with u > 1/2, 
instead, the Lagrangian dispersion is considered as super-diffusion. A well-known example 
is the Richardson's law for the particle pair separation in 3D turbulence, for which u = 3/2 
(Richardson 1926; Monin and Yanglom 1975). 

On the other hand, in the limit of infinitesimal trajectory perturbations, much smaller 
than the characteristic lengths of the system, the evolution of the particle pair separation is 
characterized by the Lyapunov exponent (Lichtenberg and Lieberman, 1982), such that 

A = lim lim - In (11.4) 

If A > the growth is exponential and the motion is said chaotic. Chaos is a familiar 
manifestation of non linear dynamics, leading to strong stirring of trajectories (Ottino, 1989). 
The process, for example, of repeated filamention around the polar vortex is basically due to 
Hamiltonian chaos (Legras and Dritschel, 1993). For finite perturbations within a smooth 
fiow, the properties of exponential separation are observed for a finite time. 



B. Finite-Scale Lyapunov Exponent 

The idea of FSLE (Aurell et al. 1997; Artale et al. 1997), was formerly introduced 
in the framework of the dynamical systems theory, in order to characterize the growth of 
non- infinitesimal perturbations (i.e. the distance between trajectories). If 6 is the scale of 
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the perturbation, and {t{6)) is the mean time that S takes to grow a factor r > 1, then the 
FSLE is defined as 

A(5) = 7^1nr (II.5) 

The average operation is assumed to be performed over a large ensemble of realizations. For 
factors r not much larger than 1, X{6) does not depend sensitively on r. If r = 2 then (t(5)) 
is also called doubling time. Operatively, + 1 scales are chosen to sample the spatial range 
of perturbation, 6o < 6i < ... < 6j\f, and a growth factor r is defined such that 6i = r ■ 
for i = 1,N. Let Imin and Imax be the smallest and the largest characteristic length of 
the system, respectively. If 6o <^ Imin then the FSLE characterizes the doubling time of 
infinitesimal perturbation. In the opposite side of the range, if 6n- ^ Imax then the FSLE 
follows the scaling law of diffusion A(5) ~ for 6 6^, as can be deduced by noticing 
that the mean square particle distance must grow linearly in time, see pi.2|) . In general, if 
the mean square size of a tracer concentration follows the (i?^) ~ t"^" law, the FSLE scales as 
X{6) ~ 5-1/^ As we have seen before, for standard diffusion u = 1/2 while for Richardson's 
super-diffusion u = 3/2. The main interest of FSLE is to select processes occurring at a 
fixed scale. We stress that definition pi.5j) differs substantially from 



defined in terms of the mean square relative displacement, because of the different averaging 
procedures in the two cases: {R"^) is computed at fixed time while t{6) is computed at fixed 
scale. As a result, a physical situation which is well characterized in terms of FSLE, either 
for scaling properties or the existence of transport barriers, may be less easily characterized 
by studying the time growth of trajectory separation (Boffetta et al., 2000b; Joseph and 
Legras, 2002). One reason is that (i?^(t)) depends on contribution from different regimes, 
as seen, for example in 3D turbulence where a dramatic dependence of R^(t) upon i?^(0) is 
observed, even at very large Reynolds number (Fung et al, 1992) 

In cases where advection is strongly anisotropic, e.g. in presence of a structure like 
the stratospheric jet stream, it may be useful to define the FSLE in terms of meridional 
(cross-stream) displacement only: 

Xmer (5^"''^ ) = , ,J lu r (II. 7) 

where S^"^^^") is the latitude distance (or meridian arc) between two points. 



Informations about the relative dispersion properties are also extracted by another fixed- 
scale statistics, the Finite-Scale Relative Velocity (FSRV), named by analogy with FSLE, 
that is defined as 



S, that is for |x(^) — x^^^l = 5. The FSRV can be regarded as the 2nd order structure 
function of the Lagrangian velocity difference and provides a complementary analysis to the 
FSLE diagnostics. In particular, in the regime of Richardson's super-diffusion, the expected 
behavior for the FSRV is ^2(6) ~ 5^/^ 

We report in the next section the results of our analysis. 

III. ANALYSIS OF THE EOLE LAGRANGIAN DATA 

After a preliminary data check, the number of balloons selected for the analysis has been 
reduced to 382. This has been obtained by discarding ambiguous ident numbers (some ident 
numbers have been used twice during the campaign), discarding trajectories that cross the 
equator and short tracks of less than 10 points. 

Successive points along a balloon trajectory were mostly recorded at a time interval of 
10~^ day (2.4 hours), but the overall distribution of the raw data does not cover uniformly 
the time axis. Hence, each of the coordinates (longitude and latitude) of every balloon 
trajectory has been interpolated in time by a cubic spline scheme, with a sampling rate of 
25 points per day. Because of possible data impurities, each Lagrangian velocity value is 
monitored at every time step (0.04 day) and data segment with abnormally fast motions are 
discarded. 

As pointed out by ML, a way to measure the dispersion between balloons is waiting 
for one of them to get close to another one, at a distance less than a threshold Sq, and 
then observing the evolution of their relative distance in time. This procedure is repeated 
for each balloon trajectory until the whole set of pairs is analyzed. The dataset includes 
original pairs of balloons that were launched within a short time interval and chance pairs 
of balloons meeting suddenly after a number of days. For the largest values of the threshold 





where 
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used in this study, the number of chance pairs largely exceeds the number of original pairs. In 
this way, global properties of the Lagrangian transport are extracted from the contributions 
of balloon pairs randomly distributed all over the domain. The number of balloon pairs and 
its evolution as the separation crosses the scales defined above is described by Table 1. 

In Figure Q four global relative dispersion curves are plotted, referring to four different 
initial thresholds Sq = 25 km, 50 km, 100 km and 200 km. The statistical samples vary 
roughly in proportional way with Sq. Relative dispersion depends sensitively, as expected, 
on the initial conditions; the four curves meet together for separation larger than about 
2500 km and saturation begins for separation larger than 4500 km, leaving room only for 
a short standard diffusive regime between these two separations and over a time duration 
of less than 10 days. The eddy diffusion coefficient, De, estimated by fitting the linear law 
4:DEt, results in De — 2.9 10^ m^s~^, a value compatible with what was found by ML. The 
pre-diffusive regime is not very clear, we can say that the behavior of the balloon separation 
looks like a power law with exponent (changing in time) between 3 and 1. 

We report in Figure El the mean logarithmic growth of the balloon relative separation 
over all pairs selected by the threshold 25 km. At very short times (< 1 day) the slope 
corresponds to an exponential growth rate with e-folding time ~ 0.4 day that we consider as 
a rough estimate of the LLE. At later times, the slope gradually decreases as the separation 
growth tends to a power law regime. In the same figure we also show the mean logarithmic 
growth of the inter-balloon distance computed for two 4-element clusters (that we label as 
'A' and 'B'), launched with a time interval of 3 days between them. A linear behavior 
(exponential growth) for both clusters is observed for short intervals; we observe as the 
exponential regime lasts longer for the 'A' cluster (~ 3 days) than for the 'B' cluster (~ 1 
days). This illustrates the fact that the duration of a dispersion regime, here the chaotic 
one, may exhibit large fluctuations generally due to different meteorological conditions. As a 
result, average time-dependent quantities, like (-R^(t)), sample different regimes at once and 
are poor diagnostics of dispersion properties. Incidentally, in ML the behavior of the relative 
dispersion between 100 and 1000 km is fitted by means of an exponential with characteristic 
e-folding time ~ 2.7 days (see ML, figure 8), which is compatible with the growth rate of 
Figure 121 between the two horizontal lines marking the range 100-1000 km, if one wants to 
fit it with an exponential curve branch. 

Figure El shows the global FSLE relatively to the same four initial thresholds used for the 
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relative dispersion and setting the amplification ratio r to -\/2- The main result of this study 
is that up to about 1000 km there is evidence of Richardson's super-diffusion, compatible 
with a k~^/^ spectrum, displayed by the behavior \{5) = ad^"^^^. The best fit is obtained 
for the initial thresholds 100 and 200 km which encompass a much larger number of pairs 
than smaller thresholds (see Table 1). The quantity is physically related to the mean 
relative kinetic energy growth rate (for unit mass) between balloons moving apart. Standard 
diffusion is approached at scales larger than 2000 km. The value of the eddy diffusion 
coefficient is estimated by fitting the FSLE in the diffusive range with (41nr)DE(5^^, as 
shown in Boffetta et al. (2000a) by means of a dimensional argument. We find that this 
value is De — 10'' m^s~^. Notice that the initial threshold does not affect very much the 
general behavior, except for obvious changes in the statistical samples. 

Figure m shows global (mainly zonal) and meridional (Amer, see pi.7|l ) FSLE of the balloon 
pairs with initial threshold 100 km. We find that the dispersion is basically isotropic up 
to scales of about 500 km, which is in rough agreement with the results of Morel and 
Larcheveque (they give a value three times larger but their analysis, see their Fig. 7, does 
not display a well-defined cut-off). At scales larger than 500 km, the two components of the 
FSLE decouple and the meridional dispersion rate follows the standard diffusion law ~ 
with a meridional eddy diffusion coefficient De ^ 10^ m^s~^. 

In order to compute the FSRV, the relative velocity between balloons is approximated 
by the finite difference formula {\R{t + At)| — |i?(t)|)/At, where \R{t)\ is the absolute value 
of the great circle arc between two balloons at time t and At = 0.04 day is the time interval 
between two successive points along a trajectory. The properties of the Lagrangian relative 
velocity are shown in Figure El The FSRV confirms the results obtained with the FSLE: 
between 100 and 1000 km the behavior is ~ 5^^^, corresponding to the Richardson's law; 
asymptotic saturation sets in beyond this range (fully uncorrelated velocities). 

IV. DISCUSSION AND CONCLUSIONS 

We have revisited the dispersion properties of the EOLE Lagrangian dataset using a 
new approach, using Finite-Scale Lyapunov Exponent, that is better suited to analyze scale 
dependent properties than standard tools that were used, e.g., by ML in a previous study of 
the same dataset. We were motivated by the fact that ML found results supporting a 
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inertial range between 100 and 1000 km, whereas more recent studies based on aircraft data 
found a k^^/^ behavior in the same range of scales. 

Our main result of our improved analysis is that the EOLE dataset supports a k~^f^ 
behavior in the range 100-1000 km as shown by the scaling properties of FSLE in this range 
indicating Richardson's superdiffusion. At distances smaller than 100 km, our results suggest 
an exponential separation with an e-folding time of about one day, in rough agreement with 
ML. At scales larger 1000 km, the dispersion tends to a standard diffusion before saturating 
at the planetary scale. Since the large-scale flow is dominated by the meandering zonal 
circulation, estimated diffusion coefficient is 10 times larger for total dispersion (10^ m^s~^) 
than for meridional dispersion (10^ m^s~^). 

Our result is compatible with an inverse 2D energy cascade in the range 100-1000 km 
or with the recently proposed alternative of a direct energy cascade (Lindborg and Cho, 
2000). Our study of the EOLE experiment has shown that this still unparalleled dataset 
of Lagrangian trajectories in the atmosphere is in agreement with results obtained from 
aircraft data. The challenge is now to compare these trajectories with the global wind fields 
produced by the recent reanalysis by operational weather centers. 
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FIGURE CAPTIONS 



Figure 1. Mean square balloon separation. The four curves refers to 4 different initial 
thresholds: 25 km (a), 50 km (b), 100 km (c) and 200 km (d). All the curves but for 25 
km have been shifted in time in order to collapse together for {E?) > 10'' km^. The eddy 
diffusion coefficient corresponding to the indicated slope is De — 2.9 10^ m^s~^. Units: time 
in days, in km^. 

Figure 2. Mean logarithmic growth of the balloon separation (+). The initial separation 
is < 25 km. The two clusters 'A' (*) and 'B' (x) have 4 balloons each and were launched 
with a 3-day lag in November 1971. The e-folding time of the exponential growth is ~ 0.4 
day. The two horizontal lines mark the range 100-1000 km (natural logarithm units). The 
straight line ML is the result found by Morel and Larchveque (1974) in their fig. 8. Units: 
time in days, R in km. 

Figure 3. FSLE of the balloon pairs, four curves with the same four initial thresholds 
as in Figure n The eddy diffusion coefficient is De — 10'' m^s~^. The quantity gives 
the order of magnitude of the relative kinetic energy growth rate (for unit mass) between 
balloons in the Richardson's regime ~ Units: 5 in km, A in day~^. 

Figure 4. FSLE of the balloon pairs, describing total (— ) and meridional (x) dispersion, 
with initial threshold 100 km. The meridional FSLE is Xmer defined in pi.7|) . The meridional 
eddy diffusion coefficient is De — 1-5 10^ m^s~^. Units: 6 in km, A in day~^. 

Figure 5. FSRV of the balloon pairs for initial threshold 50 km. The reference velocity 
is uq =100 km/hour. The slope 2/3 corresponds to the Richardson's law. Units: 6 in km. 
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TABLE I: Number of balloon pairs analyzed for each scale 
during the computation of the FSLE. The first column is the 
order of the scale (5„ = r'^So, with r = \/2; the second, third, 
fourth and fifth column refer to the initial thresholds 6o=25 
km, 50 km, 100 km and 200 km, respectively. 



n 


25 km 


50 km 


100 km 


200 km 





495 


1037 


2025 


3979 


1 


344 


782 


1670 


3649 


2 


391 


867 


1806 


3699 


3 


414 


895 


1855 


3764 


4 


440 


951 


1877 


3722 
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442 


955 


1892 


3687 
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456 


950 


1857 


3563 
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442 


944 


1829 


3471 
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448 


928 


1749 


3327 
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440 
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1703 


3131 


10 
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865 


1617 


2648 


11 
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845 


1511 
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794 


1277 
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389 


756 
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368 


642 
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346 








16 


290 
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Fig. 2: 
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Fig. 5: 
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